The Zak phase and the existence of edge states in graphene 
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We develop a method to predict the existence of edge states in graphene ribbons for a large class of 
boundaries. This approach is based on the bulk-edge correspondence between the quantized value of 
the Zak phase Z(k\\), which is a Berry phase across an appropriately chosen one dimensional Brillouin 
zone, and the existence of a localized state of momentum k\\ at the boundary of the ribbon. This 
bulk-edge correspondence is rigorously demonstrated for a one dimensional toy model as well as for 
graphene ribbons with zigzag edges. The range of k\\ for which edge states exist in a graphene ribbon 
is then calculated for arbitrary orientations of the edges. Finally, we show that the introduction 
of an anisotropy leads to a topological transition in terms of the Zak phase, which modifies the 
localization properties at the edges. Our approach gives a new geometrical understanding of edge 
states, it confirms and generalizes the results of several previous works. 
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I. INTRODUCTION 



The physics of edge states in two-dimensional (2D) sys- 
C^ems has emerged as a very challenging problem in solid 
^^tate physics. A beautiful illustration occurs in graphene, 
^ monolayer crystal of carbon,^ where the existence of such 
^tates was predicted^'^ in 1996 and confirmed experimentally 
pjater in graphene^'^ and graphene-like structures.^ This re- 
.^_Siarkable feature led to a strong activity during the last few 
^ears. For instance, edge states were predicted to give rise 
go a novel type of magnetic ordering^ and may lead to the 
'r ealization of novel spintronic devices. ^'^ 
y In a broader context, edge states are also known to play 
^n important role in quantum Hall systems^' and topo- 
QjDgical insulators. ^^'^^ Because of their chiral character, the 
"""^ge states in quantum Hall systems are robust against all 
]<^ind of disorder or interactions, while those in topological 
j-jjisulators survive scattering that preserves the time reversal 
QQymmetry. This robustness against weak perturbations can 
Q3 )e understood from a correspondence between the number of 
^^dge states and the value of a bulk topological number which 
"^s basically the Berry curvature integrated over the space of 
Qparameters, that is the Brillouin zone of the 2D system^^'^^. 

Edge states in graphene differ from those mentioned above, 
^nrthe most important distinction being that their existence 
^T^epends on the boundary conditions fixed by the shape of 
"^he edge^'^. Then, two questions naturally arise. The first 
V^ne is related to the bulk-edge correspondence if the 
r^calization of a state at the edge depends on the boundaries, 
^ it possible to relate its existence to a topological quantity 
defined within the bulk ? The second question is simply 
whether we can predict the existence of edge states for an 
arbitrary type of edge (see e.g. some examples on Fig. 1). 

These two stimulating questions have already led to many 
works. In particular, Ryu and Hatsugai showed that edge 
states in 2D systems with chiral symmetry can be related to 
a bulk topological number defined in a reduced (ID) space 
of parameters. More precisely, these authors were able to 
characterize edge states for three simple types of boundaries, 
namely the zigzag, armchair and bearded edges (a bearded 
edge is a zigzag edge with dangling bonds, also called Klein 
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FIG. 1. Examples of edges of different ribbons studied in this 
paper. They are characterized by a translation vector T(m,n), 
see section UIB. The values of the couples (m,n) are specified 
in the figure. The edges of the ribbon (1,-2) shows dangling 
bounds since mn < 0. Different ribbons can be obtained from 
the same vector T(m,n) in two different ways. For instance, the 
lattice vectors C connecting the left and right edges of the two 
ribbons defined by T(3, 3) are different. An other possibility is to 
draw a different "unit-cell" pattern for the same vector T, as for 
the two ribbons defined by T(2, 5). 
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defects^^'^^). Of course, even without disorder, there is an 
infinite number of different edge geometries in a honeycomb 
lattice, and several recent theoretical works addressed the 
existence of edge states for more sophisticated boundaries. 
Several of these works consist in tight-binding calculations 
of ribbons band structures. ^'^'^^^^ A general study including 
many various shapes of edges was also achieved by Akhmerov 
and Beenakker.^^ This work, performed within the contin- 
uous (Dirac) approximation, notably provides an analytical 
formula for the density of edge states. 

The aim of this work is to precise the bulk-edge correspon- 
dence in graphene and to address a new method to predict 
the existence of edge states for a large class of edges. We 
show that it is possible to define in an unambiguous way a 
topological phase from the bulk Hamiltonian of graphene, 
which properly takes into account the shape of the edges. 
For a ID system, this phase, called Zak phase, is nothing 
but the integration of the Berry connection over the first 
Brillouin zone^^ 



dq {Uq\dqUq) 



(1) 



where the \uq) are the Bloch wave functions. In a 2D system, 
one difficulty is to define properly the path over which the 
integration is performed, and to relate this quantity to the 
nature of the edge. In particular, for a translation invariant 
system in one direction, the Zak phase depends on the crystal 
momentum /cy associated with this direction. Here, we show 
that Z{k\\)/7r gives the number of states localized at the edge 
of the system. 

The outline of the paper is as follows. As the Zak phase 
is defined as a one-dimensional integral of the Berry con- 
nection, we first focus in Sec. II on a one-dimensional toy 
model for a chain of dimers. For this case, we give a simple 
demonstration of the bulk-edge correspondence between the 
Zak phase and the existence of edge states. Next, in Sec. Ill, 
we turn to graphene where a similar demonstration is per- 
formed for zigzag edges. By a formal analogy with the chain 
of dimers, we assume the generalization of this bulk-edge cor- 
respondence for the other boundaries. Then, we show how 
to define in an unambiguous way the Zak phase in graphene 
according to the nature of the edge and propose a very sim- 
ple graphical method to evaluate it. This information then 
directly gives us the range of k\\ for which edge states exist. 
The analytical results are in perfect agreement with those 
of Akhmerov and Beenakker^^ and reproduce many numer- 
ical works. Finally, in Sec IV we extend our approach by 
considering non-equal hopping parameters in the honeycomb 
lattice. We discuss the existence of edge states in this case 
and explain recent numerical calculations in terms of the Zak 
phase. 



II. THE ZAK PHASE AND THE EDGE STATES IN 
THE CHAIN OF DIMERS 

To illustrate the relation between the Zak phase and 
boundary states, we consider in this section a simple model, 
a one-dimensional chain of dimers A-B as shown in Fig. 2. 
The two atoms of the dimers are coupled by a hopping pa- 
rameter and the chain is obtained by coupling periodically 
the dimers with a hopping parameter t. The lattice spacing 
(between two consecutive identical atoms) is qq. 
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FIG. 2. Chain of dimers A-B. The chain starts with an atom A 
and ends with an atom B. t and t' are the hopping parameters 
and ao is the lattice spacing. The unit cell m is represented by a 
rectangle. 



The aim of this section is to show in a simple way that 
the Zak phase Z is governed by the ratio t^t, and that the 
topological transition {Z = 0) {Z = 7t) corresponds to the 
emergence of edge states in the finite system. 



A. The bulk Hamiltonian 



The Hamiltonian of the dimer chain is given by 



H = 



M 



(2) 



where aj^ (resp. creates a particle on the site A (resp. 
B) of the m}^ dimer. 

For periodic boundary conditions, we can use the Bloch 
theorem and rewrite H as 



(-f<n<f), (3) 



with *I = (Vi,fc,V'kfc)=^"'/'E^!f=ie*"°™'(«l,^l) and 

(4) 



-H^(fc) = -t 



p{k) 
p*{k) 



where p{k) = t'/t + e Introducing a = {(Tx,(Ty) the 

vector of Pauh matrices, ^^{k) can be expressed in the form 



n'^{k) = -tg{k)-a 



(5) 



with g{k) = (Re p, — Im p) = {t^t -\- cos ka^^ sink a^). Diag- 
onalizing 1-L^(k) we obtain the eigenvalues 



ek^± = ±t Uk)\=±t \p{k)\ 



= ±^t'^ ^t''^ ^2U' cos{kao) , 
and writing p{k) = \p{k)\e~'^^^^\ we have 



g{k) = \p{k)\ 



cos (j){k) \ 
sin (j){k) J 



with the phase (j){k) given by 



cot (j){k) — 



t'/t 
sin kao 



- cot kao . 



(6) 



(7) 



(8) 



The winding of the vector g(/c) as k varies across the Bril- 
louin zone is shown in Fig. 3 for two values of f /t. When 
t^t > 1, the curve g{k) does not enclose the origin and 
\(j){k)\ < 7r/2 for all k. When t^t < 1, the loop encloses the 
origin and the phase (j){k) can take any value. This topologi- 
cal behavior of the phase (j){k) is furthermore closely related 
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FIG. 3. Two trajectories of the vector g(/c) with different topolo- 
gies when k runs across the Brillouin zone. When t^t < 1 
{t' /t > 1) the curve g{k) does (not) wind around the origin. 



to the value of the Zak phase. Indeed, the eigenvectors of 
1-L^{k) are of the form 



1 / g-i0(/c) 

71 V ±1 



and the definition Eq. (1) of the Zak phase gives 

dcj) A4> 



dk- 



dk 



(9) 



(10) 



where is the variation of ^(/c) when k varies across the 
full Brillouin zone. The Zak phase Z is tt times the winding 
number of the curve g{k) around the origin, and is therefore 
zero if this curve does not enclose the origin and tt if it does. 
Thus 



Z = when t' /t > 1 
Z = TT when t' /t <1 . 



(11) 



Eq. (11) shows that tuning the ratio t' /t induces a topological 
transition (at t' /t = l) characterized by the Zak phase Z = 

Oi > Z = 7T. 



B. Open boundary conditions 

i. The missing bulk states and the Zak phase 

Consider now a finite chain of M dimers with open bound- 
ary conditions. We impose that the wave function vanishes 
at the nearest site outside the chain, that is on the B site at 
m = and the A site at m = (M + 1). Most of the eigen- 
vectors can be constructed as linear combinations |v/c,±) of 
the bulk eigenfunctions with opposite momentum \uk,±) and 
\u-k,±)- Writing \uk,±) as 



l^/c,±) 



1 

2M 



M 



m=l 



±1 



(|m,A),|m,5)) , 
(12) 

where \m^A/B) denotes the orbital A/B in the cell m, we 
obtain from the condition (0,5|v/c,±) = that |v/c,±) = 
^ [|^/c,±) - |^-/c,±)]. Thus, using (/)(-/c) = -0(/c), we have 

I^^'±)-7^Z^( ±sin(teom) ) -(l^.^).!^.^)) • 

V ^ \ / 

(13) 

Finally, the boundary condition ((M + l)^A\vk^±) = im- 
poses the quantization condition 



A;(M + l)ao - ^(fc) = A^TT, K.= 1, •••M, 



(14) 



which has to be solved in the range < k < ir/ao (the 
wave functions corresponding to k = and k = ir/ao are 
identically zero). The function (j){k) is plotted in Fig. 4, and 
the solutions of Eq. (14) correspond to the intersection of 
(l){k) with the M lines f^{k) = (M + l)kao - nir. 




kao 



FIG. 4. Variation (/)(/c) for t' /t = 1.1 (bottom curve), t' /t = 0.8 
(upper thick curve), and t' /t = 0.95 (thick dashed curve). The 
straight lines are lines of equation /^(/c) = (M + l)/cao — a^tt. 
Here, we have chosen M = 10, for which {f /t)c = 0.9091. The 
extreme values of n are indicated on the figure. In the latter case 
t' /t — 0.95 > it' /t)c^ so that there are M bulk states although 
the Zak phase Z — i^. 

From this figure, we see that the Zak phase controls the 
number of bulk states and therefore the existence of edge 
states. Indeed, when t' > t, ^(7r/ao) = and Eq. (14) has 
M solutions. When t' < t, (j){7r/ao) = tt and in this case 
there are M or M — 1 solutions depending on the value of 
t^t. By comparing the slopes of the curves (f){k) and /m(^), 
we immediately obtain that the number of bulk states \vk^±) 
depends on the critical value of the ratio t' /t 



1 



M + 1 ' 



(15) 



Including the factor 2 associated with negative and positive 
energies for each solution of the quantization condition (14), 
the number TVbuik of bulk states is: 



TVbuik = 2M 
TVbuik = 2(M 



1) 



when t'/t > {t'/t)^ 
when t'/t < {t'/t)^ 



(16) 



As we show below, the missing states are edge states localized 
at the ends of the chain. In the large M limit, the number 
of bulk states is related to the value of ^(7r/ao) since, in this 
limit, there are 2M bulk solutions when (j){7r/a) = and 
2(M — 1) bulk solutions when (piji/a) = tt. This criterion 
can be rewritten in terms of the Zak phase: since (j){k) is an 
odd function of /c, we have simply 



Z = 



dk 







dk—- — (j){7r/ao) = or tt (17) 
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As a result, by comparing Eqs. (11), (15) and (16), we con- 
clude that in the large M limit, 



^bulk 
^bulk 



2M when Z = 

2(M - 1) when Z = n . 

2. The edqe states 



(18) 



We now briefly describe the structure of the edge state 
for t' /t < {f /t)c' We search for a solution k of the form 
k = 7r/ao + ^A, where ^ = 1/A is the localization length 
of the edge state. The solution vanishing on the B site at 
m = is of the form 



M 



where 



(-1) 



m+l 



i\m,A),\m,B)) 



(19) 



(^j sinh Xaom — sinh Xao{m — 1)^ /|p(A)| | 
± sinh Xaom J 



(20) 



and has an energy 

eA,± = +t'2 _ 2tt' coshAao = ±t \p{X)\ . (21) 

The inverse localization length A is flxed by the condition 
that the wave function on the site A at m = (M+l) vanishes, 
leading to 



sinh A(M + l)ao = t sinh XMao 
Inserting (22) into (21), we flnd 

sinh Xao 



±t 



sinhA(M + l)ao ' 



(22) 



(23) 



and the components of the edge states wave functions can 
be rewritten as 



sinh Xao{M + 1 — m) 
± sinh Xaom 



(24) 



which satisfy properly the boundary conditions. Far from 
the transition, that is when the localization length ^ is much 
smaller than the size of the chain Mao, Eq. (22) reads 



1 



exp(— Aao) . 



(25) 



This implies cosh(Aao) {t^ + t''^)/2tt' , and thus eA,± ^ 0. 
More precisely, from Eq. (23), the energy vanishes as eA,± — 
exp(— AMao). The dependence X{t' /t) for a chain of M = 10 
dimers is shown in Fig. 5, together with the approximate 
expression Eq. (25). The energy spectrum of the same chain 
is displayed in Fig. 6. 



C. Remarks on the chiral symmetry 

We flnish this section with a few brief remarks concern- 
ing the Zak phase and symmetries. As flrst stressed by Ryu 
and Hatsugai,^^ the fact that the edge states have zero en- 
ergy is associated with the existence of a chiral symmetry of 
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FIG. 5. Inverse localization length A as a function of the param- 
eter t'/t, for a chain of M = 10 dimers. The full curve is the 
solution of Eq. (22). As expected, it diverges as the ratio t/t' 
reaches the critical value given by Eq. (15). The dashed line cor- 
responds to the approximation (25) valid far from the transition 
and corresponding to e = 0. 




FIG. 6. Energy levels as a function of the parameter t'/t, for a 
chain of M = 10 dimers. There is an edge state with an energy 
close to when t' /t < 1 — 1/(1 + M) . The dashed lines correspond 
to the gap ±(t^ — t) in the limit M ^ oo. The thick (dark, 
purple online) curve is the energy of the lowest energy state which 
becomes an edge state (light, green online) when (t' /t)c < 1- 



the bulk Hamiltonian.^^ From an algebraic point of view, a 
chiral symmetry is represented by an operator C which an- 
ticommutes with the bulk Hamiltonian and which satisfles 
= 1. As the Bloch Hamiltonian of the chain of dimers 
can be written as a linear combination of the Pauli matrices 
ax and ay, it is clear that C = az fulfllls these properties. 
Now, we note that the Zak phase measures the solid an- 
gle drawn by the pseudo-spinor in the Bloch sphere when k 
spans the Brillouin zone. As long as the Hamiltonian does 
not have a component proportional to a^, the pseudo-spinor 
is forced to evolve on the equator of the Bloch sphere, and 
the Zak phase is necessarily a multiple of tt. We notice that 
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breaking the inversion symmetry of the chain, for instance 
by adding a staggered potential, would add a term propor- 
tional to (Jz in the Hamiltonian and therefore would break 
the chiral symmetry. As a consequence, the Zak phase is not 
expected to be a multiple of tt anymore in this case and the 
energy of the edge states can be different from zero. 

Adding a term proportional to the identity trivially breaks 
the chiral symmetry and shifts the energy while leaving the 
pseudo-spinor on the equator so that the Zak phase is still 
quantized as a multiple of tt. A configuration where this 
simple mechanism leads to an interesting physics is obtained 
by coupling chains of dimers (assumed oriented along the x- 
direction) by a hopping parameter t" along the y direction, 
as displayed in Fig. 7. For periodic boundary conditions, the 

9 1 ^9' 1 9 1 ^9' ? 9' f 



6 • ^i' • ^6' • 6 • 6 • 
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FIG. 7. Chains of dimers coupled in the y direction by a hopping 
parameter t" . 

bulk Hamiltonian of the coupled chains reads 

Hf, = -t" cos(/c^6o)l + , (26) 

with 6o the distance between two chains. We see that the 
term proportional to 1 involves a dispersion along the y di- 
rection and does not change the properties of the Zak phase 
which are encoded in . Therefore, two distinct topologi- 
cal phases (^ = and Z = it) arise when the criteria given 
in Eq. (16) is satisfied. This means that when = tt, a 
ribbon of finite width in the x direction and invariant by 
translation in the y direction supports dispersive edge states 
along its edges, as shown in Fig 8. 
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FIG. 8. Energy levels of a ribbon built from chains of M = 
20 dimers. The chains are coupled one to each other in the y 
direction with a hopping parameter t" = 0.1, and the ribbon is 
invariant by translation along this direction. We took t = 1, and 
(left) t' = 1,5, (right) t' = 0,5. In the first case the Zak phase 
is >2^ = and there is no edge state. In the other case, the Zak 
phase IS Z — n and dispersive edge states have emerged in the 
gap. 



This example shows that the Zak phase may help charac- 
terizing the edge states even in the absence of a chiral sym- 
metry. In addition, it also illustrates that the Zak phase, 
which until now we have defined for a one-dimensional sys- 
tem, may provide informations about the edge states in sys- 
tems of higher dimension. In the following of the paper, 
we investigate the more complex case of a monolayer of 
graphene. 



III. THE ZAK PHASE AND THE EDGE STATES 
OF GRAPHENE RIBBONS 

Turning now to graphene, we demonstrate in this sec- 
tion that the Zak phase, introduced in section II for a one- 
dimensional system, has a natural generalization for a large 
class of two-dimensional graphene ribbons. In the particular 
case of zigzag edges, following the same lines as in the ID 
case, we prove that it is possible to relate the Zak phase to 
the existence of edge states. We then consider the case of 
ribbons with arbitrary orientations, and show that, comput- 
ing the Zak phase, we can predict the presence or absence of 
edge states according to the nature of the edge. 



A. The bulk Hamiltonian 

We describe the electronic spectrum of graphene by a 
tight-binding model on the triangular Bravais lattice with 
two atoms {A and B) per unit cell, as illustrated in Fig. 9. 
The parameters ti, t2 and represent the three hopping in- 
tegrals between nearest neighbors and for now we consider 
the isotropic case ti = t2 = ts (the anisotropic case is treated 
in Sec. IV). The vectors ai and a2 form a basis of the Bra- 
vais lattice, and we note (a^,a2) (with a* • slj = 27rJ^j) the 
associated basis of the reciprocal space. For periodic bound- 
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FIG. 9. (Left) Unit cell (dimer) A-B of the graphene sheet with 
the basis vectors of the Bravais lattice ai and a2 and the hopping 
parameters ti, t2 and ts. (Right) Example of an edge obtained 
by translating the dimer along ai , then twice along a2 . This edge 
is characterized by the periodicity vector T = ai + 2a2. 

ary conditions in both x and y directions, Bloch theorem 
leads to the bulk Hamiltonian 

i/^(k) = -t3|p(k)|(^,,% j (27) 
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in the basis of the two sub-lattices A and with 

p{k) = 1 + ^ exp {-ik • ai) + ^ exp (-ik • a2) 

ts ts (28) 

= |p(k)|e-'^('^) . 

The eigenenergies of the bulk Hamiltonian H^(k.) consist 
in two bands, given by e±(k) = ilit3|p(k)|, and the corre- 
sponding eigenvectors have the form, 



1 



-j0(k) 



(29) 



The positive and negative bands touch linearly at two in- 
equivalent points D and (the Dirac points) which, in the 
isotropic case considered in this section, are located at the 
corners K and of the first Brillouin zone. 

The Hamiltonian (27) can be written in the same form 
(k) = —ts g(k) • a as for the chain of dimers (5), the 
difference being that k is now a two-dimensional vector and 
that the k dependence of g(k) is of course different. As in our 
ID toy model, topological properties of the wave function as 
well as some characterization of the edge states are expected 
to be encoded in the loops drawn by g(k) as k varies across 
the Brillouin zone. This connection was actually already 
suggested by Ryu and Hatsugai^^ in the broader context of 
systems with chiral symmetry (and more recently by Mong 
and Shivamoggi in a general study of Dirac Hamiltonians^^). 
Their approach, which basically consists in the graphical 
evaluation of the Zak phase in the same way as in Fig. 3, 
was applied for three different regular types of edges (zigzag, 
armchair and bearded edges). Comparing with numerical 
calculations, they showed that the Zak phase could correctly 
predict the existence of edge states in these cases. However, 
their approach relied on the construction of a specific bulk 
Hamiltonian (i.e. a vector g(k)) for each type of edge con- 
sidered, and is therefore not convenient to consider arbitrary 
boundary conditions. Here, we keep the same bulk Hamil- 
tonian but we associate to each type of edge a specific 2D 
Brillouin zone. This allows us to predict the existence of 
edge states for different ribbon geometries, and therefore to 
address a significantly larger class of ribbons. 



B. Edges of graphene ribbons 

Turning now to graphene ribbons, we assume that both 
edges of the ribbon are parallel (i.e. that one edge can be 
deduced from the other by translation of a lattice vector C) 
and constructed as illustrated in Fig. 9, i.e. by connect- 
ing dimers A-B (or unit cells) of fixed orientation (vertical 
orientation in Fig. 9). 

More precisely, considering two positive or negative inte- 
gers m and n, an edge is built as (|m| + |n|) translations of 
the dimer, \m\ of which along ai (— ai if m negative) and 
|n| of which along a2 (— a2 if n negative), in an arbitrary 
order, and by repeating the pattern obtained in this way. 
Therefore the edges are invariant under the translation vec- 
tor T = mai + na2 which characterizes the type of edge. 
Noting the angle between T and the y axis, this angle is 
related to (m, n) through 



When m and n have the same sign, the class of edges con- 
structed in this way exactly corresponds to the minimal 
boundary conditions of Akhmerov and Beenakker.^^ In the 
other case, the edges exhibit dangling bonds. 

Fig. 1 gives some examples of such edges. It is easy to see 
that the vectors T(m, n) and T(n, m) describe the same kind 
of edge. We also notice that the same vector T can describe 
edges with different shapes, see the example of T = (2,5) in 
Fig. 1. 

We stress that choosing the "unit-cell" dimer A-B with a 
different orientation (i.e. rotated from zb27r/3 with respect 
to the vertical one used in Fig. 9) leads to a different set 
of boundaries. These latter are of course just deduced from 
the former by a ±27r/3 rotation, and thus considering only 
the vertical unit-cell dimer, as we shall do in the following, 
does not restrict the type of edge to be studied. We insist 
however that if one wanted to consider another orientation 
of the unit-cell dimer for the edges, it would be essential 
to modify accordingly the dimer orientation in the definition 
of the bulk Hamiltonian (which basically fixes the zero of 
the phase 0(k) in Eqs. (27) or (29)). This is a necessary 
condition to derive a bulk-edge correspondence for the edge 
states in terms of the Zak phase. 



C. The Zak phase in graphene ribbons 

For one-dimensional models such as the one considered in 
section II, the Zak phase is defined as the integral of the 
Berry connection across the Brillouin zone. To generalize 
this notion to two-dimensional systems such as graphene, 
this integration should be taken on a cut of a 2D Brillouin 
zone in a direction transverse to the ribbon orientation. More 
precisely, as the ribbon is assumed to be invariant under 
translation by the vector T, Bloch theorem guarantees that 
/cy, the component of the crystal momentum parallel to T, 
is a good quantum number. We expect therefore the Zak 
phase 2{k\\) to be a function of /cy, and to correspond to an 
integration of the Berry connection across the 2D Brillouin 
zone along a perpendicular direction k±. 

Let us assume for now that m and n are coprime inte- 
gers (we will return below to the case where they are not). 
We choose the Brillouin zone from which the Zak phase will 
be computed as the one generated by two orthogonal vec- 
tors of the reciprocal space, Fy and F^, obtained as fol- 
lows. The first of these vector Fy = 27rT/|Tp is parallel to 
the direction T of the ribbon and merely defines the one- 
dimensional Brillouin zone of the ribbon. The second vector 
is taken perpendicular to T, and its norm is fixed by 
the constraint that Fy x F^ = x ^2 (where, as mentioned 
above, a^ and a.2 are the reciprocal lattice vectors defined by 
a^ • a^ = 27rSij). This leads to 



Fy(m,n) = 



(n + 2m)a^ + (m + 2n)a^ 



2(n2+m2 
r^(m, n)= na-l — ma.2 . 



nm) 



(31) 
(32) 



tan^ 



1 n 



y/3n-\-m 



(30) 



For the sake of completeness we remind briefly in appendix A 
why (Fy , F^) constructed in this way actually defines a Bril- 
louin zone when (m, n) are coprime integers, and only in this 
case. 

For an arbitrary edge characterized by the vector T(m, n), 
we introduce the unit vectors ey = Fy/lFyl and e_L = 
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r^/|r^|, and write the momentum as k = k\\e\\ + k±e±. 
The Zak phase Z{k\\) can thus be defined as: 

Z(k\\) =i j) dk^{u\^^±\dk^u\^,±) , (33) 

which, by using the expression of the Bloch function (29), 
simply reads 

Z{k\\) = ]^jdk^dkM\^) . (34) 

Note that as is a vector of the reciprocal lattice, the 
integration can indeed be seen as taken on a closed path. 

Let us consider now the situation where m and n are not 
coprime, in which case (Fy , F^) obtained from Eqs. (31) and 
(32) do not form a basis of the reciprocal lattice (see ap- 
pendix A). Writing m = Irh and n = In with / integer 
and in and h coprime and following exactly the same line 
of argument as above, we can construct a basis (r||,r^) 
of the reciprocal lattice corresponding to (m,n), and define 
the Zak phase ^(rh,n) (^||) accordingly from (34). An ex- 
ample of such a construction will be shown in Sec. HIE 
(see the case (m, n) = (2,0) of Fig. 13). However, as 
T(m,n) = T(m,n)//, we have now Fy = /Fy. As a con- 
sequence, to a given value k\\ of the quantum number of the 

ribbon correspond / values (k^^\ • • • , Brillouin 

zone (^||-^^ = k\\ + j|f|||), and therefore / Zak phases. The 
prescription we shall use for non-coprime integers (m, n) is 
therefore that for a given value k\\ of the ribbon quantum 
number, the Zak phase is defined as 

I 

^(m,n) (^11 ) = ^Z{7h,n){k\\^) • (35) 

i=i 

This prescription merely corresponds to a folding of the Bril- 
louin zone. 



D. Zigzag boundary conditions 

Before we address general orientations for the graphene 
ribbon, let us first consider in details the simple case of 
zigzag boundary conditions. In that case, and as illustrated 
in Fig. 10, the two edges of the ribbon, (Si) and (S2), are 
constructed as the vertical dimer A-B translated periodi- 
cally with the vector T = ai (see Fig. 9). We note C the 
vector which connects the two edges and C = C + 2a2 — ai 
the vector connecting B'^ and lines of empty sites 

nearest neighbors to the edges Bi and 62- Noting ^k± = 
(^J^^, ^kdr)^ Bloch state of momentum k, we have 

^k,±(r) = (r|^k±)a^(^%^ j . (36) 

Following the approach used for the one-dimensional chain of 
dimers of section JIB, we would like to construct eigenstates 
of the ribbon as a linear combination ^ of two Bloch states 
^kzt and ^k'± at the same energy. From Fig. 10, we see 
that the boundary conditions read 



y 




FIG. 10. Schematic representation of a zigzag ribbon. The zigzag 
edge is obtained by translating the dimer A-B (represented by a 
rectangle) with the periodicity T = ai. The two edges Bi and 
B2 are represented by thick lines. The first vacant sites outside 
of the ribbon where we impose the wave function to vanish are 
represented by circles. C is the vector of the Bravais lattice which 
connects two sites on both edges, and connects the first vacant 
site on one side of the ribbon to the first vacant site on the other 
side. 

with G Z. 

The boundary condition Eq. (38) implies that the wave 
functions ^ are combinations of the form ^ = — ^k' 
with k • T = k' • T (here T = ai). Since we need ek = 
ek', we have to consider momentum pairs (k, k') such that 
(k + k') • a2 = k • ai = k' • ai. For a given value of k\\ , it can 
be checked that this is satisfied if 

k = k|| + kx , (39) 

= k|| - kx , (40) 

where we have introduced ky = /cyey and k^ = k^e^. Note 
that k^ = or k^ = F^/2 correspond to ^k — ^k' = 0- 
The boundary condition Eq. (38) then imposes 

(k-kO • C - (0(k) - 0(kO) = 2ktt (41) 

(for integer n)^ or in term of the phase ^(/cy, k^) = ^[^(ky + 

k^)-(/>(ky-k^)] 

k^-C'-0(A:y,A:^) = /^7r, (42) 

which is the strict equivalent of the quantization condition 
(14) obtained for the ID chain. We note furthermore that 
(Fx • C) = 2tt{M + 1), with M the number of dimers in the 
transverse direction of the ribbon, and thus half the number 
of bands. Using the same arguments as in section II B, we 
see that the number of edge states depends on the value of 
(/)(/cy , Fx/2), and more precisely that the quantity 




gives the number of pairs of edge states (with opposite ener- 
gies) for a given k\\. 

E. General orientation 



*'^(C' + i/T) = [on B'^] , (37) 
*^(i^T)=0 [onBi], (38) 



We consider now graphene ribbons characterized by an 
arbitrary translation vector T(m, n). Based on similar argu- 
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merits as in the previous section we conjecture that the re- 
lation between the Zak phase and the number of edge states 
holds in this general case. Then we use this relation to pre- 
dict the existence of edge states in graphene ribbons of gen- 
eral orientation. We also check that for every case for which 
edge states has been computed (numerically or otherwise) 
their appearance is correctly predicted by the Zak phase. 

In Eq. (34), the phase ^(k) should be understood as a 
multivalued function. The single-valued function (^(k) corre- 
sponding to the restriction of 0(k) to the interval [— tt, +7r] is 
displayed on Fig. 11. It shows lines of discontinuity connect- 
ing pairs of Dirac points. The location of the discontinuities 
(i.e. which Dirac points are connected by them) depends on 
the choice made for the unit cell dimer A-B (oriented along 
the y axis in this paper, see Fig. 9). Therefore the phase 
(/)(k) is not invariant by a rotation of an angle ±27r/3. 



















Sr 




\ 

\ 


1 






1 1 1 1 r 







FIG. 11. Density plot of the phase 0(k). The discontinuities of 
the phase are shown by horizontal white lines connecting pairs 
of Dirac points. These singularities separate the values cj) — +7r 
(light region) to the value = — tt (dark region). Consequently, 
the only paths that contribute to a non-vanishing Zak phase 
are those that cross these singularities. The black lines, thick 
lines and dashed lines represent the iso-0 lines respectively for 
= 0, = -7r/2 and +7r/2. 

Let us start again by considering ribbons of period T(m, n) 
with coprime {m^n). The discontinuities of 0(k) are ex- 
tremely convenient to determine the Zak phase since paths 
that lead to a non-vanishing Zak phase necessarily cross a 
discontinuity line. Actually, the Zak phase >2'(m,n) (^||) is 
given by the number of discontinuities d{k\\) intersected by 
the path V^rn,n)iJ^\\) = [k||,k|| + T^] along which Z(^,^)(A:||) 
is computed, that is 



(44) 



In other words, d{k\\) is just the number of pairs of edge 
states (of opposite energies) for a given /cy . 

These considerations make it possible to compute graph- 
ically the Zak phase in a rather straightforward way. For a 



given choice of edge characterized by the vector T(m,n) = 
mai + na2, we first represent the vector T±{m^n) = na^ — 
ma.2 (Eq. 32). This is done in Fig. 12, where we plot the 
vectors T±{m^n) from the left extremity of the discontinu- 
ity (0,0) to the left extremity of the discontinuity (n, — m). 
Next, in Fig. 13, we translate perpendicularly T± until the 
left extremity of an another discontinuity is reached (dashed 
line), which gives T\\{m^n). The rectangle defined by T± 
and Fy is the Brillouin zone we want to associate with the 
ribbon. For a given value of the momentum /cy, 2{k\\) is 
then deduced from Eq. (44) by simply counting the num- 
ber of intersections of the segment [k||,k|| + T±] with the 
discontinuity lines of ^(k). 




FIG. 12. On top of the density plot of the phase 0(k), we have 
represented several values of the vector T±{m,n) = na.1 — ma.2. 
The left extremity of a discontinuity of 0(k) is taken as the origin. 
We have explicitly plotted the vector r^(l, 0) as an example. 

If m and n are not coprime, as for instance in Fig. 13(d), 
r||(m,n) and r^(m, n) cannot be obtained in this way 
as they do not define a Brillouin zone. Writing however 
(m, n) = /(m,n) with (m, n) coprime, the Brillouin zone 
corresponding to (m,n) can be constructed as above, and 
one has simply T\\{m^n) = T\\{m^h)/l and T±{m^n) = 
lTj_{rh^h), which basically amounts to a folding of Fy by 
a factor /. For a given value of /cy, the Zak phase in then 
given by Eq. (35). 



F. Range of existence and density of edge states 

In this section we derive the range A/cy for which edge 
states exist, as well as the density of edge states for arbitrary 
boundary conditions thanks to the bulk-edge correspondence 
in terms of the Zak phase. 

In the last section, we showed that the number of pairs of 
edge states for a given value of k\\ is given by the number 
d{k\\) of crossings between the path 7^(m,n) (^||) and the dis- 
continuities of 0(k). Since a Brillouin zone always contains 
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FIG. 13. Surfaces |r|| A r^| associated with different ribbon 
vectors T(m,n). T± is obtained as T±{m,n) = na.1 — ma.2 = 
(n, — m); see text for the construction of r||. (a) T±(m = l,n = 
2) = (2,-1), (b) ri(-2,3) = (3,2), and (c) Fx (0,1) = (1,0). 
In these three cases, m and n are coprime, and the surfaces 
|r|| AF^I represented by a shaded rectangles are Brillouin zones, 
(d) T(0, 2) : the surface obtained by Fn and F^ is not a Brillouin 
zone because m and n are not coprime. The corresponding Bril- 
louin zone (shaded rectangle) is given by F^(0, 1) = F^(0,2)/2 
and F||(1,0) = 2F||(2,0). 



exactly one line of discontinuities, the total range 

rW|T| 



d{k\\)dk\\ 



(45) 



over which the ribbon exhibits edge states, is obtained by 
projecting the line of discontinuities onto the k\\ axis, as il- 
lustrated in Fig. 14. This leads to 



47r 

Ak\\ = ^Isin^l 

Otto 



(46) 



where is the angle between the direction T of the ribbon 
and the vertical axis y of the dimers (see Fig. 9 and Eq. (30)). 
Therefore, there is no edge state for edges parallel to the 
armchair edge {6 = 0) and Ak\\ is maximum and equal to 
47r/3ao for bearded edges. ^^'^^ Then, from Eq. (30) we have 
the relation between the range A/cy of existence of edge states 
in graphene and the integers (m, n) characterizing the edge: 



Ak\\ (m, n) 



27r 



n — m 



3ao ^/\ 



- nm 



(47) 



We notice that Ak\\ (m, n) = Ak\\ (Im^ In) which means that 
T(m,n) and lT{m^n) support the same number of local- 
ized states. A relevant quantity to study is the ratio IZ = 




FIG. 14. Projection of the discontinuity of (/)(k) onto the k\\ axis. 
This projection gives the range AA;|| for the existence of edge 
states. 



AA:||/|r||| = A/c|| |T|/(27r) which gives the relative range of 
the ID Brillouin zone where edge states exist. We find 



7^ = 



(48) 



We now comment the formulas (47) and (48). First, we 
apply these formulas to several types of edges listed in table 
I. The three general types of edges mentioned in this table 





(m, n) 


A/cii 


n 


armchair 


(1,1) 








zigzag 


(1,0) 


27r/3ao 


1/3 


bearded zigzag 


(1,-1) 


47r/3ao 


2/3 


bearded armchair 


(2,-1) 


27T 1 A/3ao 


1 




(1,3) 


47r/3yi3ao 


2/3 


general types 


(1,5) 


27r/3V3lao 


4/3 




(1,8) 


147r/15V3ao 


7/3 



TABLE I. Applications of the formulas (47) and (48) for some 
edge geometries studied in previous works. 

represent a zigzag profile extending over respectively three, 
five and eight unit cells followed by an armchair defect (see 
for example Fig. 1 for the case (1,5)). These results are in 
good agreement with the size of the edge states energy bands 
obtained by numerical tight-binding calculations in different 
previous works. ^'^'^^'^^^^ We also notice that in the limit 
\n\ ^ |m| (or the other way), we recover A/cy 2i\ jZa^^ 
which is the expected result for the zigzag edge. 

The formulas (47) and (48) directly lead to the impor- 
tant result that edge states exist for most types of ribbons 
with a periodic pattern, which is in agreement with a previ- 
ous analytical approach within the Dirac framework. More 
precisely, we find that there is no edge state if and only if 
n = m. As already mentioned, this class of ribbons includes 
all the vectors T parallel to the y axis of the dimers A — B 
(see for instance the case (3,3) displayed in Fig. 1. Similar 
ribbons have been synthesized for the first time recently^^). 
The well-known particular armchair case corresponds to the 
smallest |T| that obeys this condition. 

The case IZ > 1 also deserves some attention. It im- 
plies that two or more pairs of localized states may corre- 
spond to the same k\\. This situation happens when, as 
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illustrated in the example of Fig. 13(b), several disconti- 
nuities are intersected by the path 7^(m,n) (^||)- Such a situ- 
ation is automatically achieved for m and n coprime when 
6{m,n) > 6^(4,1) = arctan(^/3/5). Moreover, by construc- 
tion, the projection of the line of discontinuity of (j) spans 71 
times the Brillouin zone, which implies that the number of 
localized states cannot differ by more than one unit for any 
two k\\ (see Fig. 15). Therefore, for m and n coprime, there 



E 













A"A 



k|| 



k|| 



k|| 



FIG. 15. Schematic band structures of graphene ribbons that 
exhibit edge states. Edge states are represented as zero-energy 
flat bands as expected in the large width limit within the tight- 
binding model (assuming the chiral symmetry is preserved). The 
number of pairs of edge-states is indicated on the figure, (a) 
corresponds to the case IZ < 1, (b) corresponds to 7^ > 1. (c) Our 
analysis based on the Zak phase predicts that this configuration 
is not allowed. 

is either [7^] or [7^] + 1 edge states for each k\\ {[x] is the floor 
function). In the particular case where 1Z is an integer there 
are exactly IZ edge states for each /cy . 

Finally, we can define the quantity p = 7l/\T\ which cor- 
responds to the "density of edge states per unit length" in- 
troduced by Akhmerov and Beenakker. We get 



27T 



3ao 



nm 



(49) 



which was first obtained in Ref. [23] for minimal boundary 
conditions by a different method. 



IV. EMERGENCE AND DESTRUCTION OF EDGE 
STATES: A TOPOLOGICAL APPROACH 

In this section, we generalize the formula (47) for non- 
equal hopping parameters ti ^ t2 7^ ^3, and establish a cri- 
terion for the existence of edge states for an anisotropic hon- 
eycomb lattice. We show that the manipulation of these pa- 
rameters leads to a topological transition described in terms 
of Zak phase that affects the range A/cy of existence of edge 
states. We stress that breaking the isotropy of the hopping 
parameters preserves the chiral symmetry, and therefore the 
topological character of the Zak phase. As a consequence 
the analysis developed in the previous section generalizes 
straightforwardly to the case considered below. 



A. Effect of an anisotropy on the existence of edge 
states 

Several previous works dealing with the tight-binding 
model in the honeycomb lattice showed that the Dirac points 
move when modifying the ratio of the parameters ti/tj?^^'^ 



As the lines of discontinuities of (/>(k) connect pairs of Dirac 
points, the modification of the ratio ti/tj changes the Zak 
phase and therefore leads to new ranges A/cy of existence of 
edge states. This is clearly shown in Fig. 16. To determine 




FIG. 16. Density plot of the phase 0(k) for t2 = 1.5ti = l.Sts- 
The discontinuities of the phase are shown by curved white lines 
that pair the Dirac points. The black lines, thick lines and dashed 
lines represent the iso-0 lines respectively for = 0, = — 7r/2 
and +7r/2. 

these new ranges, we have first to specify the position of the 
Dirac points for any ti. Up to a vector of the reciprocal 
lattice, they are given by 



with 



D 

d2 =^e 



27r 

TT — 

27r * 



arccos 



arccos 



27r 
27r ^ 



(50) 



/2 



2t2t3 
+ tl- 



2hh 



(51) 



where Re(x) takes the real part of x. The cartesian coordi- 
nates od the Dirac points D^'^ = {Dx \ Dy^) are given by 



D 



(^{di + ^2 - 27r) /ao, {di - ^2) /VSao^ 
(^{-di - d2^ 27r) /ao, {d2 - di) /V^a^ . 



(52) 



Then, since the Dirac points D^'^ are not located anymore 
at the corner K(') of the Brillouin zone, the range A/cy of 
existence of edge states is modified as: 

AA:|| = |D-D'||sin((9 + /3)| (53) 

where P is the angle between the line [D'D] and the kx axis 
(see Fig. 17), which is then given by 

27r/ao -I^; 



cos/3 =2- 



sin/3 =2- 



|D-D^ 
D-D1 



(54) 
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FIG. 17. Projection of the discontinuity line of 0(k) onto the k\\ 
axis. The asymmetry of the hopping induces a modification of 
the discontinuity locations as compared with Fig. 14. 

Using Eqs. (30) and (54), the range for the edge states 
reads: 



AA:ii = 



(n - m) (l^ - I^;) + (n + m) D^^ 



■ (55) 



Next, by using the expression of the positions D^^ and Dy 
given in (52) one finds: 



Ah 



as well as 



2 \nd2 — mdi \ 



\nd2 — mdi\ 



(56) 



(57) 



where di and d2 are given by Eqs. (51). These results 
give a criteria for the existence of edge states that links the 
anisotropy of the hopping parameters encoded in di and d2 
with the nature of the edge characterized by {m^n). In the 
isotropic case, we have di = d2 = 7r/3, and we recover the re- 
sult (47) discussed in the previous section. The formula (56) 
means that an edge state exists in graphene-like structures 
if \ndi — md2\ 7^ 0. 

An interesting consequence is that, for a given type of rib- 
bon, edge states can emerge or collapse when an anisotropy 
is applied. For instance, edge states can emerge for armchair- 
like boundary conditions (m = n) when either ti/ts ^ 1 or 
^2/^3 7^ 1 (see Refs. [25] and [33] for the case n = m = 1). 
This is clearly displayed in Fig. 18. In the same way, edge 
states can collapse by manipulating the asymmetry of the 
hopping in such a way that \ndi — md2\ = 0. 



B. Merging of Dirac points and edge states 

As one increases the anisotropy of the system, for instance 
by modifying one of the ratio ti / tj , one may eventually reach 
a point where ^1=^2+ ^3 (or the equivalent up to a cyclic 
permutation of the indices). From Eq. (51), this condition 
implies that di and ^2 take the values or tt, which, from 
Eq. (52), corresponds to a merging of the Dirac points. By 
increasing further the anisotropy, a gap opens at the merg- 
ing point. This merging is a topological transition since the 
Berry phases ±7r associated to the two Dirac points annihi- 
late at the transition. ^^'^^ 



1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 
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FIG. 18. (Left) Density plot of the phase 0(k) for ti = 1, t2 = 1.5 
and ts = 1 represented in the Brillouin zone corresponding to 
armchair edges. The two horizontal lines delimit the region where 
2{k\\) = TT. (Right) Band structure of an armchair ribbon with 
the same hopping parameters. In this case, k\\ = ky. The two 
vertical lines delimit the same range as in the left panel, which is 
such that edge states at zero energy, clearly separated from the 
bulk bands, have emerged. 



This topological transition is a bulk property, and is thus 
independent of the orientation of the anisotropy. On the 
other hand, the existence of edge states at the merging tran- 
sition and beyond depends on the orientation of the ribbon 
with respect to this anisotropy, and this information is still 
contained in the Zak phase. Indeed the anisotropy of the 
hopping parameters controls the size and the location of the 
lines of discontinuities of 0(k). It is therefore essential to dis- 
tinguish which pair of Dirac points merges. There are three 
possibilities : 

• ^3 > ti + ^2 

In this configuration, the Dirac points that merge are 
the endpoints of the same discontinuity line. Then, the 
discontinuities disappear at the merging {ti + ^2 = ^3) 
and the Zak phase vanishes. This leads to 



7^ = 



(58) 



whatever m and n (see top panel of Fig. 19), imply- 
ing that edge states never exist in this case. As an 
illustration, we plot in Fig. 19 the band structure of a 
zigzag ribbon exactly at the merging point: the zero 
energy edge-states^ of the isotropic case have collapsed 
because the Zak phase is zero for all /cy . 

The situation is totally different when the two merging 
Dirac points are attached to two distinct discontinuity lines, 
which occurs when ti — t2 = ±^3. In this case, the two 
discontinuity-lines themselves merge, implying that IZ is an 
integer, as we discuss now. 

• ti>t2^t3 

In this case we have d2 = and di = which leads to 



7^ = Iml 



(59) 



for all n. For instance, there is no edge state for zigzag 
edge dit = +7r/6 (see center panels of Fig. 19 ), but 
there is one over the whole ribbon Brillouin zone for 
armchair edges and for zigzag edges at ^ = — 7r/6. We 
notice that the existence of edge states now only de- 
pends on the orientation of the edge given by but not 
on k\\ anymore. 
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• ^2 > ti + t3 

In this case di = and d2 = tt and we find 

n = \n\ (60) 

whatever the value of m. Of course this case is equiv- 
alent to the one discussed previously with the substi- 
tution m ^ n that is ^ —0. The phase (^(k) and 
the corresponding band structure for zigzag ribbons at 
6 = +7r/6 are represented in Fig. 19 (bottom panels). 

In summary, we note that in all three cases, at the merg- 
ing, the Zak phase becomes independent of /cy and remains 
unchanged in the gapped phase. This implies that at the 
merging transition and beyond, 1Z is necessarily an integer. 
This is obvious from Figs. 19 : the merging of the Dirac 
points implies either a disparition of the discontinuity lines 
or their transformation into an infinite line. 




-10 -5 5 10 



kxCiQ k||ao 

FIG. 19. Density plot of the phases 0(k) at the merging of the 
Dirac points, and band structure for zigzag ribbon with = 7r/6 
(that is T(m = 0,n = 1) ^ T±{m = 0,n = 1) = (1,0)). (Top) 
= 1, ^2 = 1 and ts — 2, (center) ti = 2, t2 = 1 and ts = 1, 
(bottom) ti = 1, t2 = 2 and ts = 1 . The scale is given in units 
of 1/ao. 



compute the Zak phase in order to predict the number of 
edge states. The approach we have developed consists in 
constructing the appropriate 2D Brillouin zone associated 
with the vector T(m, n) which defines the edge. The Zak 
phase Z{k\\) giving the number of edge states for each /cy is 
then directly obtained by integrating the Berry connection 
along a path fixed by T(m,n) and k\\ in this 2D Brillouin 
zone. 

We stress that this bulk-edge correspondence is, beyond 
the ID chain of dimers case, only rigorously proven here for 
zigzag edges. It is therefore so far a conjecture for the class 
of edges we have defined. This conjecture is however sup- 
ported by the fact that it reproduces all the known previous 
results obtained (numerically or otherwise) in the literature 
for various specific types of edges. ^'^'^^'^^^^ 

In practice, the value of the integral defining the Zak phase 
is easily obtained graphically. Our approach therefore does 
not require any sophisticated formalism or calculation, and 
gives an elegant understanding of the origin of the edge- 
dependent edge states in terms of a topological bulk quan- 
tity. In particular, it provides a simple understanding of the 
appearance and disappearance of edge states by manipulat- 
ing the anisotropy of the tight-binding hopping parameters. 
Such a manipulation may be induced in graphene by apply- 
ing an uniaxial stress or bending of the sheet, or in pho- 
tonic crystals which mimic the same physics, by changing 
the distance between the confining mirrors.^ 

We finish with a few comments concerning the connection 
between Zak phase and edge states. 

First, this bulk-edge correspondence differ from the ones in 
quantum Hall systems^^^^ or Z2 topological insulators^^'^^ 
since here the existence of edge states precisely depends on 
the orientation of the edge. This difference with the usual 
bulk topological numbers originates from the fact that the 
Zak phase is a ID (rather than 2D) integral of the Berry 
connection. 

Second, we stress that within our approach, the vector 
T(m, n) defining the periodicity of the ribbon entirely deter- 
mine the Zak phase. As many different shapes may corre- 
spond to the same vector T(m, n), the Zak phase and there- 
fore the number of edge states are expected to be indepen- 
dent of the variation of the edge geometries as long as they 
correspond to the same T(m,n).^^ 

Finally, our description of edge states in terms of the Zak 
phase is a priori not restricted to graphene but should be 
in principle also applicable to other 2D systems like d-wave 
superconductors,^^ square lattice with half a quantum flux 
per unit cell^^ or bi-layer graphene^^ for instance. 
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V. CONCLUSION 

In this paper, we have investigated the correspondence 
between the Zak phase and the existence of edge states for 
arbitrarily oriented graphene ribbons with a large class of 
edge shapes. We have proposed a definite prescription to 



Appendix A: Construction of the Brillouin zone 

In this appendix, we give a brief reminder of the reason 
why the vectors Fy and T± defined by Eqs. (31)- (32) actually 
generate a Brillouin zone when n and m are coprime integers. 

This latter condition indeed implies that one can find two 
integers (m^, n') such that mn' — nm' = 1, in which case the 
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couple of vector (T, N), with N = m'ai +n'a2, form a basis 
of the Bravais lattice. The choice of (m', n'), and thus of N, 
is not unique, but this is irrelevant for our purpose. 

From (T, N), one deduce a basis (r^v, T^) of the reciprocal 
lattice, 

Tn = n'al -m'sil ^ (Al) 
Fx = nai - , (A2) 

which is such that T±_ _L T. A Brillouin zone can thus be ob- 



tained from the parallelogram generated by (Fat, F^). More 
generally however, any vector F such that (F — F^v) || F^ is 
such that the parallelogram generated by (F,F^) is a Bril- 
louin zone. A natural choice is to take for F the vector Fy 
which is parallel to T (and thus orthogonal to F^). Since 
T . Fiv = 27r(mn' - nm') = 27r, one has |F|| | = 27r/|T| which 
is nothing but the size of the (ID) Brillouin zone of the rib- 
bon. 
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